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ABSTRACT 

We investigate the stellar content of the starburst dwarf galaxy IC10 using accurate and deep optical 
data collected with the Advanced Camera for Surveys and with the Wide Field Planetary Camera 2 
on board the Hubble Space Telescope. The comparison between theory and observations indicates a 
clear change in age distribution when moving from the center toward the external regions. Moreover, 
empirical calibrators and evolutionary predictions suggest the presence of a spread in heavy element 
abundance of the order of one-half dex. The comparison between old and intermediate-age core He- 
burning models with a well defined overdensity in the color-magnitude diagram indicates the presence 
of both intermediate-age, red clump stars and of old, red horizontal branch stars. 
Subject headings: galaxies: individual (IC10) — galaxies: stellar content — Local Group — stars: 
evolution 



1. INTRODUCTION 

Dwarf galaxies are ubiquitous stellar systems outnum- 
bering giant systems in the Local Group (LG, Mateo 
1998), in the Local Volume (d < 10 Mpc, Vaduvescu & 
McCall 2008), and in the nearby Universe (Popesso et 
al. 2006; Milne et al. 2007). Recent evidence indicates 
that the Local Group includes at least 62 dwarfs and 
among them 26 ± 5% are dwarf irregulars (dlrrs, Grebel 
2003; McConnachie et al. 2008). However, we still lack 
firm criteria discriminating between dlrrs and Blue Com- 
pact Galaxies (BCDs). According to Thuan (1985) and 
to van den Bergh (2000) the BCDs are dlrrs that are ex- 
periencing a significant burst of star formation. On the 
other hand, Richer & McCall (1995) found that the metal 
abundance of BCDs is more similar to dwarf sphcroidals 
than to dlrrs, and Papaderos et al. (1996) pointed out 
the lack of an evolutionary link among BCDs, dlrrs and 
dwarf ellipticals (dEs). This key issue is far from being 
settled, and indeed in a recent detailed photometric and 
spectroscopic investigation Vaduvescu & McCall (2008) 
suggested that BCDs, dlrrs and dEs define the same fun- 
damental plane. 

Dwarf irregulars also play a key role in constraining the 
impact that structural parameters and intrinsic proper- 
ties have on the evolution (initial mass function, star for- 
mation history) of complex systems (Massey et al. 2007). 

1 Based on observations collected with the ACS and the WFPC2 
on board of the HST. 

2 Univ. Roma ToV, via della Ricerca Scientifica 1, 00133 Rome, 
Italy; Nicoletta.Sanna@roma2.infn.it 

3 INAF-OAR, via Frascati 33, Monte Porzio Catone, Rome, 
Italy 

4 DAO-HIA, NRC, 5071 West Saanich Road, Victoria, BC V9E 
2E7, Canada 

5 INAF-OACTe, via M. Maggini, 64100 Teramo, Italy 

6 IAC, Calle Via Lactea, E38200 La Laguna, Tenerife, Spain 

7 STScI, 3700 San Martin Drive, Baltimore, MD, 21218, USA 

8 Univ. Pisa, Largo B. Pontecorvo 2, 56127 Pisa, Italy 

9 INFN, Sez. Pisa, via E. Fermi 2, 56127 Pisa, Italy 

10 INAF-OAT, via G.B. Tiepolo 11, 40131 Trieste, Italy 

11 ESO, Karl-Schwarzschild-Str. 2, 85748 Garching bei 
Munchen, Germany 

12 CTIAO-NOAO, Casilla 603, La Serena, Chile 



Moreover, they are fundamental laboratories for investi- 
gating the evolution of massive stars in systems that are 
undergoing significant bursts of star formation (Crowther 
2007). 

Among the nearby dlrrs IC10 is an interesting sys- 
tem, since it is one of the most massive (log M/M =8.49, 
Vaduvescu et al. 2007), and the comparison between the 
H a and the _B-band luminosity indicates that it is expe- 
riencing a starburst phase (Hunter et al. 1993). More- 
over, it has been suggested that IC10 harbors a large 
number of young Wolf-Rayet stars (Massey & Holmes 
2002) and intermediate-age carbon stars (Demers et al. 
2004). However, we still lack detailed knowledge of 
the stellar content of IC10. In particular, Vacca et al. 
(2007), using deep optical and near-infrared data, pro- 
posed that the isochrone fit to IC10 — at fixed distance 
and metal content — would require different reddening 
values for Main Sequence (MS) and Red Giant Branch 
(RGB) stars. 

In a previous investigation (Sanna et al. 2008) we 
provided a new estimate of the distance modulus 
(/i=24.60±0.15 mag) based on a new calibration of 
the tip of the RGB, and of the reddening (E(F555VF- 
F 1 814H / )=1.16±0.06 mag) based on empirical calibra- 
tors. Here we address the galaxy's stellar content. 

2. RESULTS AND DISCUSSION 

The photometric catalog we adopt is based on archival 
optical images from the Advanced Camera for Surveys 
(ACS) and the Wide Field Planetary Camera 2 (WFPC2; 
see top panel of Fig. 1 and Sanna et al. 2008). The 
final catalog includes ~ 720,000 stars with at least one 
measurement in each of two different bands. The ACS 
data in the F555W and F81AW bands were placed on 
the VEGAMAG system following Sirianni et al. (2005). 
The F606IF-band images collected with the ACS were 
transformed into the F^SVF-band using local standards, 
and the same approach was adopted to transform the 
^5551^ and the F8UW images collected with WFPC2 
into the corresponding ACS bands. On average the star- 
to-star precision of the above transformations is better 
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than 0.02 mag (Sanna et al. 2009, in preparation). The 
final catalog was split into two different regions: region 
C) covers the galaxy center and includes both ACS and 
WFPC2 data, while region E) lies at a radial distance 
greater than two arcminutes and only includes ACS data 
(see the blue and red polygons in the top panel of Fig. 1 
and Sanna et al. 2008). 

Spectroscopic estimates of IClO's metallicity, based 
on HII regions, indicate a metal content ([Fe/H] ~ 
-0.71 ±0.14, Lee et al. 2003) similar to the Small Magel- 
lanic Cloud (SMC, [Fe/H] « -0.7, Zaritsky et al. 1994 

also based on HII regions or [Fe/H] 0.75 ± 0.08, 

Romaniello et al. 2008, based on Ccphcids). To char- 
acterize the stellar content of IC10 we adopted differ- 
ent sets of both scaled-solar (young and intermediate 
ages) and a-enhanced (old ages) isochrones from the 
BaSTI database 13 plus a few young isochrones specif- 
ically computed for this project. In particular, we 
adopted isochrones based on evolutionary tracks ac- 
counting for mass-loss (77=0.4), neglecting both convec- 
tive overshooting during the core H-burning phases and 
atomic diffusion (Pictrinfcrni ct al. 2004, 2006). We have 
assumed a true distance modulus of /i=24.60 and a red- 
dening E(F555W-F8UW)=l.W mag. Data plotted in 
the bottom left panel of Fig. 1 show that young scaled- 
solar isochrones (red lines) at fixed metal and helium 
content ([M/ff]=-0.66, F=0.251) and ages ranging from 
6 Myr {M{Turn - Off[TO])/M Q ~ 30.5) to 100 Myr 
(M(TO)/M Q ~ 4.3) properly fit young MS stars in IC10. 
The same conclusion applies to the fit of RGB stars 
with the old a-enhanced isochronc (t<~14 Gyr). Note 
that the global metallicity of this isochrone is [M/H]=- 
0.66 with [Fe/H]=-1.0l and [a/Fe] -0.4 (Pietrinferni 
et al. 2006). The data plotted in the bottom panel 
of Fig. 1 were selected according to photometric error 
(cf814W=cf555W_ 0.1), separation (sep > 4) and sharp- 
ness (\sha\<0.3). The separation index quantifies the de- 
gree of crowding (Stetson ct al. 2003). The current value 
corresponds to stars that have required a correction of 
less than a few percent for light contributed by known 
neighbours. The sharpness index quantifies the similar- 
ity between the shape of the measured objects and of the 
Point-Spread- Function (PSF). It is the quadrature differ- 
ence between the one-sigma-half-width of the measured 
object and the one-sigma-half-width of the core of the 
PSF (Stetson & Harris 1988). The bottom right panel 
of Fig. 1 shows that the CMD of the external regions is 
deeper, since these regions are less affected by crowding. 
The comparison with the central regions indicates a sig- 
nificant change in age distribution, and indeed, massive 
MS stars almost disappear when moving toward the ex- 
ternal regions. The MS stars are well fit by isochrones 
with ages ranging from 12 Myr (M(TO)/M Q ~ 14.2) to 
350 Myr (M(TO)/M Q ~ 2.5). These findings and the 
continuous stellar distribution along the MS indicate that 
IC10 experienced ongoing star formation during w the 
last half Gyr. Fig. 1 also indicates that IClO's stel- 
lar populations show a spread in metal content. The 
width in color of RGB stars is well fit by isochrones 
with a single age (14 Gyr) and metal contents rang- 
ing from [M/H] ~ -0.96 (Y =0.248, green line) to 
[M/H] 0.35 (Y =0.256, blue line). The above com- 

13 http://www.oa-teramo.inaf.it/BASTI 



parison between theory and observations has to be con- 
sidered as a preliminary guideline. These estimates of 
age and metal content are affected by empirical (distance 
modulus, reddening, photometric zero-points) and theo- 
retical (mixing length, color-temperature relations) un- 
certainties. Firm constraints on these parameters require 
deep and accurate photometry down to the TO of the old 
population. Note that the possible presence of differen- 
tial reddening amounting to ±10% would not account for 
the observed spread in color (see the arrows in the right 
panel of Fig. 1 and Sanna et al. 2008). However, IC10 has 
an extended HI envelope, a large number of HII regions 
(Hidalgo-Gamez 2005) and molecular clouds (Leroy et al. 
2006). This means that internal spatial variations of the 
reddening are quite probable. To partially overcome this 
problem, we adopted as representative of the IC10 stellar 
content the stars located in a small external region (E, 
see the red polygon in the top panel of Fig. 1). 

To further characterize the stellar content of IC10 we 
also adopted empirical tracers. Fig. 2 shows the compar- 
ison between field E) of IC10 and the ACS photometry of 
an SMC field provided by Sabbi et al. (2007; red dots). 
Note that in this figure we plotted a number of IC10 
stars, randomly selected, similar to the number of stars 
present in the SMC field. The SMC sample was plot- 
ted by assuming for the SMC a true distance modulus of 
fi = 18.9 and a reddening of E(B -V) = 0.08 mag. The 
ages of the main stellar components in this SMC field 
range from a few tens of Myr (bright MS) to a few Gyr 
(red clump, RC). To investigate the possible presence of 
an old stellar population we also compared IC10 to the 
globular cluster NGC 6362. The V, /-band photometry 
(Stetson 2000) for this cluster was transformed into the 
ACS photometric system using the transformations by 
Sirianni et al. (2005). NGC 6362 is an almost metal- 
rich cluster ([Fe/H] = -1.04 ± 0.06, [M/H] —0.75, see 
Table 1) and its Horizontal Branch (HB) morphology 
is characterized by both red and blue stars (Brocato et 
al. 1999). Data plotted in this figure, in particular in 
the helium burning region (i.e., RC and red HB stars, 
25.8<F814T^<25.2, \.7<FbbbW - F814VF<2.1), indi- 
cate that IC10 hosts both intermediate-mass stars (RC) 
and candidate old low-mass (red HB) stars. To further 
constrain these key points Fig. 3 shows the comparison of 
the same field with three GCs with different metal con- 
tents ([M/H] — 0.47, 47 Tuc; -0.85, NGC 2808; -0.93, 
NGC 1851, see references listed in Table 1) and HB mor- 
phologies (only red HB stars, 47 Tuc; red, blue HB and 
RR Lyrae stars, NGC 1851; red, blue HB and extreme 
HB stars, NGC 2808). Note that in this figure we plot- 
ted a number of IC10 stars, randomly selected, similar to 
the number of stars present in the GC NGC 2808. Data 
in Fig. 3 further support the evidence (see Fig. 1) that 
RGs in IC10 cover at least one- half dex in metal content 
(-0A<[M/H]<-1). Moreover, a fraction of the stars lo- 
cated near F8UW^26 and 1.5<F555VF - F8UW<1.8 
appear to be candidate RR Lyrae stars. Finally, we note 
that the comparison with empirical calibrators indicates 
that both metal-intermediate and metal-rich candidate 
RGB bump stars (see Table 1) should be systematically 
fainter than our current limiting magnitude at the color 
typical of RGB bump stars (F555IF - F8UW^2. 1-2.2 
mag, see the arrows in Fig. 3). 

To investigate the properties of He-burning stars we 
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F555W-F814W [mag] F555W-F814W [mag] 

Fig. 1. — Top - The coverage of the HST data sets collected with the ACS and with the WFPC2 (black lines). The blue and red polygons, 
superimposed to IC10, mark fields C) and E), respectively. The background is a MegaCam@CFHT image of 6x7 arcmin. North is up and 
East is to the left. Bottom left - F814W, F555W-F814W CMD of IC10 central regions. Red lines display scaled-solar isochrones (BaSTI) 
at fixed chemical composition and different ages (see labeled values). The circle (M(TO)/Mq=30.5), the diamond (M (TO)/Mq=7.7) and 
the triangle (M(TO)/Mq=4.3) mark the Turn-Off (TO) of three young isochrones. Bottom right - Same as the left, but for the tCfO 
external regions. The cross (M(TO)/M =14.2), the asterisk (M(TO)/Mq=5.0) and the square (M(TO)/M =2.5) mark the TO of three 
young-intermediate age isochrones. The green and the blue lines show two old (t=14 Gyr) isochrones at different metal contents. The 
arrows in the top right corner show the shift in color and in magnitude caused by the possible occurrence of a differential reddening of 
±10%. 



adopted different sets of evolutionary models constructed 



assuming both old and intermediate-age progenitors. 
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Fig. 2.— F814W, F555VK-F814W CMD of IC10 external regions 
(black dots) compared with an SMC field (red dots) and with the 
GC NGC 6362 (green dots). The number of IC10 stars plotted 
in this figure is similar to the number of stars in the SMC field. 
See text and Table 1 for more details concerning the true distance 
moduli and the reddenings adopted for these systems. 



Fig. 3. — Same as Fig. 2, but the comparison is performed 
with three GCs: NGC 1851 (green dots), NGC 2808 (red dots) 
and 47Tuc (blue dots). The number of IC10 stars plotted in this 
figure is similar to the number of stars in NGC 2808. The true 
distance moduli and reddenings adopted for these systems are listed 
in Table 1. From top to bottom the arrows display the position of 
the RGB bumps for NGC 1851 and 47Tuc. 



The top panel of Fig. 4 shows a-enhanced Zero-Age- 
Horizontal-Branch (ZAHB, solid line) models together 
with the He-exhaustion locus (dashed, 10% of central 
He still available) at fixed metal and helium content 
([M/H]=-0.ffi, F=0.251) for an old (simU Gyr) pro- 
genitor (M pr =0.80 Mq). To cover the age range of 
IC10 stars we also adopted core He-burning (solid) and 
He-exhaustion (dashed, 10% of He left) models for a 
set of intermediate-age (160 Myr <t< 6.6 Gyr) pro- 
genitors (1.0<M pr <3.5 Mq). The evolutionary proper- 
ties of He-burning, low-to-intermediate mass stars have 
been thoroughly investigated in the literature (Sweigart, 
Greggio & Renzini 1990; Castellani et al. 2000; Pietrin- 
ferni et al. 2004,2006; Bertelli et al. 2008). Here we 
note that the ratio between He- and H-burning life- 
times is quite constant (tue/tH =0.006) when moving 
from M(HB)=0.60 to 0.80 M (old progenitor). On the 
other hand, the same ratio changes from tu e /tu =0.11 
(M(RC)=1.83M Q ) to 0.39 (M(RC)=2.18Mq) and to 
0.34 (M(RC)=2.78M Q ) U for scaled-solar, intermediate- 

14 The use of He lifetimes at He-exhaustion, i.e., no He left in 
the core, changes the quoted ratios by a few thousandths and a few 
hundredths for old and intermediate-age progenitors, respectively. 



mass progenitors (see top panel of Fig. 4). This 
stark difference is caused by the fact that when mov- 
ing from 0.8 to 2.2 M/Mq the core He-ignition takes 
place in structures that are less and less affected by 
electron degeneracy. This means that the He core 
mass at He-ignition, and in turn the luminosity dur- 
ing core He-burning, decreases from Mf Ie /MQ—0A85 
(M pr =0.80 M Q , M fsuw =-0.35 mag, t H -14 Gyr) to 
M|r e /M Q =0.467 (M pr =2.2 Mq, M F814W =-0.05 mag, 
i#=750 Myr). More massive structures are character- 
ized by a convective core during H-burning phases and a 
further increase in stellar mass causes a steady increase 
in the He core mass and in luminosity (Af pr =2.80 Mq, 
MJj e =0.370 M Q , M F814W =-1.02 mag, t H =28Q Myr). 
The above difference implies that the expected star count 
ratio between MS and He-burning structures increases 
by 1-2 orders of magnitude when moving from old to 
intermediate-mass stars. 

Therefore, we decided to perform a more detailed 
comparison between theory and observations. We se- 
lected stars in the external regions using severe crite- 
ria (<JF8iiw=o'F555W < 0.07 mag, sep > 6, \sha\<0.2). 
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Data plotted in the bottom panel of Fig. 4 show that the 
well defined ovcrdensity centered on F814W~25.72 and 
F555W - F814W~1.90 agrees quite well with the ex- 
pected position of old, red HB stars and of intermediate- 
age red clump stars. The solid lines display the ZA- 
HBs for two different metal contents ([M/H]=-0.96, 
Y =0.248; [M/H]=-0.35, Y = 0.256) and for old and 
intermediate-age progenitors. The stellar masses and the 
ages of the progenitors are quite similar to the models 
plotted in the top panel. Note that the spread in magni- 
tude and color of the He-burning region is larger than the 
typical photometric errors (see error bars). The range in 
color covered by the stellar overdensity is also system- 
atically bluer and larger than the color range covered 
by typical low-mass RGB bump stars. The RGB bump 
in a metal-rich (47Tuc) and in a metal-intermediate 
(NGC 1851) globular cluster is, indeed, fainter and red- 
der (see the arrows in Fig. 3). This comparison also sup- 
ports the hypothesis that the stars with F814W~26 and 
F555W - F814W~1.7 are candidate RR Lyrae stars. 
Note that corrections for completeness of the current 
photometry would go in the direction of increasing the 
number of candidate old HB stars. 

The photometric accuracy in the region around the 
peak does not allow us to distinguish RC from old HB 
stars. However the occurrence of warm HB stars, once 



confirmed by independent experiments, will provide a 
robust identification of the so-called Baade's red sheet, 
i.e., evidence for an old stellar population (Baade 1963) 
in a starburst galaxy. Current circumstantial evidence 
is, indeed, based on intermediate-age (RC) He-burning 
stars (Aparicio et al. 1997; Schulte-Ladbeck et al. 1998). 
Moreover, the identification of massive MS stars, old 
(HB), and intermediate age (RC) helium burning stars 
indicates that IC10 underwent several star formation 
episodes during its life. A CMD a couple of magnitudes 
deeper and with a stronger temperature sensitivity could 
provide firm constraints on whether the star formation 
activity of this interesting system has been continuous or 
sporadic. 
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Fig. 4. — Top - Color-Magnitude diagram for predicted He- 
burning structures at fixed global metallicity ([M/H]=— 0.66). The 
fainter solid and dashed lines show the a-enhanced ZAHB and 
the He-exhaustion (10% of He core still available) for an old pro- 
genitor (t H =14: Gyr, Af pr =0.80 M Q ). The brighter solid and 
dashed lines show the core He-burning and the He exhaustion for 
intermediate-age progenitors. The mass of the progenitors, the 
mass at core He-burning and the ratio between He and H lifetimes 
for selected structures are labeled and marked with black (old) 
and red (intermediate-age) symbols, respectively. Bottom - The 
red solid polygons display the 35, 60, 80 and 97% isodensity levels. 
Fainter blue and green lines show the ZAHBs for old, low-mass 
structures with different chemical compositions. The almost ver- 
tical green and blue lines display core He-burning structures with 
the same compositions, but for intermediate-mass structures. 
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TABLE 1 

Intrinsic parameters of the GCs adopted as empirical calibrators. 



ID 




E(B - V) b 


[Fe/Hf spe 


[Fe/H]% w 


[Fe/H\ d CG 


[M/H] e 


RGB bump f 


47Tuc 


13.32 ± 0.09 g 


0.04 ± 0.02 h 


-0.76 ±0.04' 


-0.71 ±0.05 


-0.78 ± 0.02 


-0.47 


13.49 ± 0.10 


NGC6362 


14.43 ± 0.05 j 


0.08 ± 0.02J 


-1.04±0.06 k 


-1.18 ±0.06 


-0.99 ± 0.03 


-0.75 


14.45 ± 0.10 


NGC2808 


15.23 ±0.10' 


0.18 ± 0.01 m 


-1.14±0.10 n 


-1.36 ±0.05 


-1.11 ± 0.03 


-0.85 


15.10 ± 0.10 


NGC1851 


15.44 ± 0.20° 


0.02 ± 0.02° 


-1.22±0.03 p 


-1.23 ±0.11 


-1.03 ±0.06 


-0.93 


15.16 ± 0.10 



a Cluster true distance modulus (mag). 
b Cluster reddening (mag). 

c High-resolution spectroscopic iron abundances. 

d Iron abundances based on Ca triplet measurements provided by Rutledge et al. (1997) in the Zinn & West (1984) and in the Carrctta 
& Gratton (1997) mctallicity scale. 
° Global mctallicity based on spectroscopic iron abundances and assuming [a/Fe] = 0.4 (Salaris et al. 1993). 

f The Johnson-Cousins magnitudes of the RGB bump provided by Di Cecco et al. (2009, in preparation) were transformed into the ACS 
F814W-band (VEGAMAG) following Sirianni et al. (2005). 

SBono et al. (2008). h Salaris et al. (2007). ^Koch ct al. (2008). ^Brocato ct al. (1999). fc Gratton (1987). 'Castellani et al. (2006). m Bedin 
et al. (2000). "Carretta (2006). "Savianc ct al. (1998). ^Yong et al. (2009). 



